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NUMERICAL  SIMULATION  OF  TRANSIENT 
HYPERVELOCITY  FLOW  IN  AN  EXPANSION  TUBE 


P.  A.  Jacobs^ 

Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  VA  23665 


ABSTRACT 

We  present  several  numerical  simulations  of  the  transient  flow  of  helium  in  an  expansion 
tube.  The  aim  of  the  exercise  is  to  provide  further  information  on  (and  hopefully  further  in¬ 
sight  into)  the  operational  problems  of  the  NASA  Langley  expansion  tube.  The  calculations 
were  performed  with  an  axisymmetric  Navier-Stokes  code  based  on  a  finite- volume  formula¬ 
tion  and  upwinding  techniques.  Although  laminar  flow  and  ideal  bursting  of  the  diaphragms 
was  assumed,  the  simulations  showed  some  of  the  important  features  seen  in  the  experi¬ 
ments.  In  particular,  the  discontinuity  in  tube  diameter  at  the  primary  diaphragm  station 
introduced  a  transverse  perturbation  to  the  expanding  driver  gas  and  this  perturbation  was 
seen  to  propagate  into  the  test  gas  under  some  flow  conditions.  The  disturbances  seen  in  the 
test  flow  can  be  characterized  as  either  “small-amplitude”  noise  possibly  introduced  during 
shock  compression  or  “large-amplitude”  noise  associated  with  the  passage  of  the  reflected 
head  of  the  unsteady  expansion. 


'This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No.  NASl-18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 
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Nomenclature,  Units 

A  :  (x,  y)-plane  cell  area 

a  :  local  speed  of  sound,  m/s 

Cp  :  coefficient  of  heat  capacity  (constant  F),  Jjkg 

Cy  :  coefficient  of  heat  capacity  (constant  volume),  J/kg 

E  :  total  energy  (internal  +  kinetic),  J/ kg 

e  :  specific  internal  energy,  J/kg 

e  :  unit  vector 

F  :  algebraic  vector  of  x-component  fluxes 

G  :  algebraic  vector  of  y-component  fluxes 

h  :  specific  enthalpy,  J/kg 

k  :  coefficient  of  thermal  conductivity 

n  :  unit  normal  vector 

P  :  pressure.  Pa 

Pr  :  Prandtl  number,  (Cpfi/k) 

Q  :  algebraic  vector  of  source  terms 
q  :  heat  flux 

R  :  gas  constant,  J/kg/K 

Re  :  Reynolds  number 

r  :  radial  coordinate,  m 

T  :  temperature,  K 

t  ;  time,  s 

U  :  algebraic  vector  of  conserved  quantities 
TI  :  Riemann  invariant,  m/s 

u  :  x-component  of  velocity,  m/s 

V  :  y-component  of  velocity,  m/s 

X  :  x-coordinate,  m 

y  :  y-coordinate  (radial),  m 

/3  :  parameter  in  Mirels’  boundary  layer  theory 

p  :  density,  kg/m^ 

7  :  ratio  of  specific  heats 

A  :  second  coefficient  of  viscosity 

T  :  shear  stress.  Pa 

p  :  coefficient  of  viscosity,  Pa.s 

n  :  cell  volume,  m® 

n'  :  volume  per  radian  for  the  axisymmetric  cell 

<  •  >  ;  cell- averaged  value 
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Subscripts 


0 

1  ...  20 


e 

i 

V 

w,  wall 


:  conditions  immediately  behind  the  shock 

:  gas  state  as  shown  on  Fig.  1 

:  conditions  just  outside  the  boundary  layer 

:  interaction  of  expansion  with  driver- gas/test-gas  interface 

:  viscous  contribution 

:  wall  value 

:  cartesian  components 
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1  Introduction 


The  expansion  tube  is  a  member  of  the  family  of  pulse-type  aerodynamic  facilities  designed 
to  provide  high  enthalpy  test  gas  for  short  durations.  It  has  a  fundamental  advantage 
over  the  reflected-shock  tunnel  because  its  operating  cycle  does  not  involve  the  stagnation 
of  the  heated  test  gas.  Hence,  it  can  provide  a  final  test  flow  with  low  ionization  and 
low  dissociation.  Initial  theoretical  studies  of  the  expansion  tube  [1]  indicated  that  such 
a  machine  could  provide  a  wide  range  of  test  flows  by  simply  altering  the  initial  filling 
pressures.  However,  operational  experience  (see  e.g.  [2,  3])  indicated  that  the  test  flow 
was  often  contaminated  by  large-amplitude  disturbances.  The  purpose  of  this  report  is  to 
present  some  numerical  simulations  of  an  entire  (but  idealized)  expansion  tube  in  order  to 
provide  some  further  information  (and  hopefully  some  insight  into)  the  mechanisms  causing 
the  disturbances  in  the  test  flow.  The  specific  facility  configuration  and  the  flow  conditions 
considered  here  are  based  on  the  experimental  perfect-gas  study  undertaken  by  Shinn  and 
Miller  [3]. 


1.1  Facility  Operation 

An  expansion  tube  is  essentially  a  single  tube  divided  into  three  sections  by  “primary”  and 
“secondary”  diaphragms  as  shown  in  lower  part  of  Fig.  1  (which  is  reproduced  from  [3]).  The 
strong  primary  diaphragm  separates  the  driver  tube  and  the  intermediate  (or  shock)  tube 
while  the  light  secondary  diaphragm  separates  the  intermediate  tube  from  the  acceleration 
tube.  The  intermediate-  and  acceleration  tubes  have  the  same  diameter.  Initially,  the  driver 
tube  contains  a  “driver  gas”  at  high  pressure  (state  4)  while  the  intermediate  tube  contains 
the  test  gas  (state  1)  and  the  acceleration  tube  contains  a  low  pressure  “acceleration  gas” 
(state  10).  The  aerodynamic  model  to  be  tested  is  located  near  the  downstream  end  of  the 
acceleration  tube. 

The  operation  of  the  facility,  starting  with  the  rupture  of  the  r.rimary  diaphragm  at  <  =  0, 
is  shown  schematically  in  the  top  part  of  Fig.  1.  The  numbering  of  the  flow  states  is  the  same 
as  that  introduced  by  Trimpi  [1].  High  pressure  driver  gas  expands  into  the  intermediate 
tube  (state  3)  and  shock-compresses  the  test  gas  (state  2).  Simultaneously  an  expansion  fan 
travels  upstream  into  the  driver  tube.  On  reaching  the  end  of  the  intermediate  tube,  the 
shock  ruptures  the  secondary  diaphragm  (ideally  without  producing  disturbances)  and  the 
leading  portion  of  the  shock- compressed  test  gas  is  allowed  to  expand  into  the  acceleration 
tube  to  reach  the  test  flow  conditions  (state  20).  The  expanded  test  gas  will  have  a  high 
velocity  (typically  6  —  8  km/s)  and  a  relatively  low  temperature  (several  hundred  degrees 
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Kelvin). 

As  shown  in  Fig.  2,  test  time  commences  after  the  passage  of  the  test-gas/acceleration- 
gas  interface  and  ideally  finishes  with  the  arrival  of  the  downstream  end  of  the  unsteady 
expansion  (and  a  corresponding  increase  in  pitot  pressure).  Not  all  of  the  shock-compressed 
test  gas  is  expanded  to  the  test  flow  state  (20)  and  so  test  times  are  typically  measured 
in  tens  to  hundreds  of  microseconds.  When  the  upstream-head  of  the  unsteady  expansion 
intersects  the  driver-gas/test-gas  interface,  a  partial  reflection  will  occur  and,  for  some  test 
conditions,  this  wave  (which  travels  along  a  u-fa  characteristic)  will  arrive  at  the  end  of  the 
acceleration  tube  before  the  unsteady  expansion  itself.  Unfortunately,  ideal  operation  of  the 
facility  is  seldom  realized  as  test  flows  are  typically  very  noisy  and  may  have  peak-to-peak 
variations  of  50%  in  pitot  pressure. 


1.2  Previous  Studies 

The  cause  of  the  test  flow  disturbances  has  been  the  focus  of  several  studies.  Shinn  and 
Miller  [3]  made  experimental  measurements  in  the  NASA  Langley  facility  using  helium  as 
both  test  and  driver  gas  in  order  to  eliminate  chemical  effects.  They  observed  that  most  of  the 
operating  conditions  resulted  in  test  flows  which  were  very  noisy  but  identified  an  optimum 
condition  (Pi  =  3.45  kPa,  Piq  =  16  Pa)  with  approximately  300  fis  test  time  before  the 
onset  of  large-amplitude  fluctuations  in  the  pitot  pressure.  A  large  dip  in  the  pitot  pressures 
was  observed  for  values  of  Pio  significantly  larger  than  20  Pa  and  this  was  attributed  to 
boundary  layer  transition.  They  also  identified  the  importance  of  the  diaphragm  dynamics 
on  the  quality  of  the  test  flow,  especially  at  low  initial  pressures,  but  they  did  not  identify 
the  mechanism  causing  the  early  disturbances  in  the  test  gas. 

More  recently,  equivalent  experiments  [4]  using  both  helium  and  argon  as  the  working  gas 
were  performed  at  the  University  of  Queensland  to  check  the  Langley  results  [3]  and  see  if  the 
useful  range  of  test  conditions  could  be  expanded.  In  conjunction  with  these  experiments, 
a  numerical  study  (based  on  quasi-one-dimensional  modelling)  was  undertaken  [5]  to  see  if 
the  test  gas  disturbances  could  be  caused  by  “blobs”  of  dense  driver  gas  being  accelerated 
through  the  expanding  test  gcis.  It  was  concluded  that,  even  if  the  blobs  existed,  they  would 
not  arrive  early  enough  to  cause  the  test-flow  disturbances.  However,  the  reflection  of  the 
unsteady  expansion  off  the  driver-gtis/test-gas  interface  could  (under  some  conditions)  arrive 
early  enough  to  introduce  the  noise.  Test  time  estimates  based  on  this  mechanism  are  given 
in  [6]  and  are  consistent  with  the  present  numerical  simulations. 

Pauli  [4,  7,  8]  has  pursued  the  possibility  of  acoustic  disturbances  being  generated  near 
the  primary  diaphragm  station  and  then  being  propagated  through  the  unsteady  expansion 
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of  the  driver  gas.  A  family  of  transverse  waves  can  exist  and,  when  propagated  through  the 
expansion,  will  undergo  a  frequency  shift  such  that  a  single  dominant  frequency  will  emerge. 
Transmission  of  these  waves  across  the  interface  and  into  the  test  gas  is  assumed  to  be 
determined  by  the  acoustic  impedance  mismatch  at  the  interface.  For  situations  where  the 
sound  speed  in  the  test  gas  is  much  larger  than  that  in  the  expanded  driver  gas,  attenuation 
of  the  waves  in  the  test  gas  should  occur.  Otherwise,  a  significant  fraction  of  the  noise  will  be 
transmitted.  This  simple  rule-of-thumb  does  not  seem  to  hold  for  the  simulations  presented 
here. 


2  The  Numerical  Method 

The  computations  reported  here  were  performed  with  a  finite-volume  upwind  code  based 
on  the  full  Navier-Stokes  equations.  The  code  is  described  briefly  here  but  further  detail  is 
available  in  [9], 

2.1  Governing  Equations 

For  an  axisymmetric  flow  (with  y  as  the  radial  coordinate),  the  finite- volume  formulation  of 
the  Navier-Stokes  equations  may  be  expressed  as 

^  l(yF  -  vF^)  dy  -  ^  -  yo.)  dx^Q'  (i) 

where  <  •  >  indicates  a  volume  average.  O'  is  the  volume  per  radian  and  5  is  a  contour  in 
the  (x,  y)-plane  around  the  volume.  The  U,  F  and  G  vectors  are 
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Note  that,  except  for  the  “y”  premultiplying  factor,  they  are  the  same  as  those  in  the  planar 
two-dimensional  situation.  The  viscous  terms  are 
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and 


^  dv  ^  ( du  dv  v\ 
r„=2;,^  +  A(^^  +  ^  +  -j  . 

(5) 

f  du  dv\ 

(6) 

CQ  1'^ 

1 

II 

CO 

1 

II 

(7) 

Treating  the  viscous  contributions  in  the  form  yr  avoids  any  difficulties  with  the  geometry 
singularity  at  2/  =  0.  The  effective  source  term  is 


0 

0 

(P  -  Tes)A/n' 

0 


,  ^  V  .( du  dv  v\ 

where  Tgg  =  2fi-  +  A  —  +  - 

y  \dx  ay  yj 


(8) 


To  augment  these  equations,  we  use  the  equation  of  state  for  a  calorically  perfect  gas 

P  =  p{l-  l)e  ,  (9) 

with  7  =  1.667  and  R  =  2077  J/kglK.  The  corresponding  specific  heats  are  =  3114  Jlkg/K 
and  Cp  =  5191  Jjkg/K.  The  first  coefficient  of  viscosity  is  evaluated  as 

H  =  5.023  X  10"^  Pa.s  ,  (10) 

and  the  second  coefficient  is  obtained  from  Stokes’  hypothesis  A  =  Also,  a  constant 

Prandtl  number  of  0.67  is  assumed. 


2.2  Numerical  Implementation 

The  flow  domain  in  the  (x,j/)-plane  is  discretized  as  a  structured  mesh  of  quadrilateral  cells 
with  flow  properties  stored  at  the  cell  centres.  At  each  time  step,  the  inviscid-flux  vectors 
(2)  are  evaluated  by  first  applying  a  MUSCL  interpolation  scheme  to  obtain  “left”  and 
“right”  states  at  the  midpoints  of  the  cell  interfaces.  A  locally  one-dimensional  (approximate) 
Riemann  solver  is  then  applied  to  obtain  the  interface  flow  properties  during  the  time  atep. 
The  spatial  derivatives  used  in  the  viscous  flux  vectors  (3)-(7)  are  obtained  at  the  cell  vertices 
by  applying  the  divergence  theorem.  The  source  term  (8)  is  evaluated  at  the  cell  centres.  The 
line  integrals  in  (1)  are  then  evaluated  using  the  midpoint  rule  and  the  solution  advanced  in 
time. 

Although  a  predictor- corrector  time-stepping  scheme  was  available  in  the  code,  the  solu¬ 
tions  shown  in  the  following  sections  have  used  Euler  time  stepping  in  order  to  reduce  the 
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required  cpu  time  for  each  simulation.  Fig.  3  shows  the  density  contours  in  the  shock  and 
driver  tubes  at  t  =  2.0  ms  for  case  2.  Except  for  the  fine-scale  features,  there  are  few  differ¬ 
ences  in  the  solutions  for  both  Euler  time-stepping  and  predictor-corrector  time-stepping  on 
the  1604  X  20  grid.  The  differences  between  these  “coarse-grid”  solutions  and  the  solution 
on  the  (fine)  2406  x  30  grid  are  also  small.  The  limit  on  the  number  of  cells  was  set  by  the 
memory  available  on  the  SUN  workstations  (64  Mbytes)  and  the  desire  to  fit  the  jobs  within 
the  8  Mword  queues  on  the  Cray-YMP  (Sabre). 


3  Simulation  Results 

Figure  4  shows  the  flow  domain  considered  in  the  simulations.  The  flow  is  considered  to  be 
axially  symmetric  with  the  flow  domain  extending  radially  from  the  axis  to  the  tube  wall 
and  axially  from  the  closed  end  of  the  driver  tube  to  the  end  of  the  acceleration  tube  where 
it  enters  the  test  section  of  the  facility.  Dimensions  are  shown  in  the  figure.  The  grid  was 
generated  in  three  pieces,  one  for  each  section  of  the  facility.  Note  that  the  number  of  cells 
tind  the  clustering  toward  the  tube  wall  was  kept  constant  along  the  length  of  the  domain. 
Also,  the  cell  size  remained  constant  in  the  axial  direction.  The  radial  distance  from  the 
outermost  cell-centres  to  the  intermediate-  and  acceleration-tube  walls  is  3.87  x  10“^  m. 

Note  that,  while  the  change  in  tube  radius  at  the  primary  diaphragm  station  is  included, 
neither  the  diaphragm  nor  the  square  transition  section  downstream  of  this  point  is  modelled. 
Instead,  the  simulations  start  with  high-pressure  reservoir  conditions  upstream  of  the  primary 
diaphragm  station  and  intermediate-tube  initial  conditions  downstream  of  the  same  station 
(and  into  the  acceleration  tube).  The  simulations  are  allowed  to  proceed  until  the  primary 
shock  reaches  the  secondary  diaphragm  station  and,  at  that  point  in  time,  the  conditions 
downstream  of  the  secondary  diaphragm  are  reset  to  acceleration-tube  initial  conditions. 
Thus,  the  bursting  of  the  diaphragms  is  considered  to  be  ideal  in  the  following  calculations. 
This  is  in  contrast  to  the  experimental  observation  that  the  flow  quality  is  very  sensitive  to 
rupture  behaviour  of  the  diaphragms  [2,  3]. 

3.1  Initial  State 

All  of  the  simulations  used  the  nominal  driver  tube  conditions  discussed  in  [3].  These  are 

P4  =  33  MPa,  U4  =0,  T4  =  330  K,  p4  =  48.15  kgjm^,  64  =  1.028  MJjkg. 

As  discussed  in  [3]  the  nominal  driver-tube  temperature  may  be  lower  than  the  actual  gas 
temperature  at  primary-diaphragm  rupture  where  a  maximum  value  of  T4  =  390  K  was 
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estimated.  The  test  gas  and  acceleration  initial  conditions  are 


T\q  —  T\  —  300  C\Q  —  Cl  —  0.934  M J f  kg ,  wio  —  'Wi  —  0, 

with  pressures  and  densities  as  specified  in  Table  1. 

For  cases  la  and  lb,  the  conditions  were  chosen  to  approximate  the  optimum  conditions 
identified  in  the  experiment  [3]  and  two  different  grids  were  used  to  check  for  sensitivity  to 
grid  resolution.  Cases  either  side  of  this  optimum  condition  were  then  considered.  Case  2  has 
a  higher  intermediate-tube  pressure  where  large-amplitude  noise  was  observed  at  an  earlier 
time  in  the  experiment.  Case  3,  with  a  low  intermediate-tube  pressure,  also  was  observed 
to  have  a  noisy  test  flow  early  in  the  test  time.  Finally,  case  4  considers  the  condition  with 
high  acceleration-tube  pressure  where  (for  the  Langley  expansion  tube)  a  large  dip  in  pitot 
pressure  was  observed  shortly  after  the  arrival  of  test  flow  [3]. 


3.2  Case  1:  Pi  =  3.45  kPa^  Piq  =  16  Pa 

This  condition  was  chosen  to  be  the  starting  point  for  the  study  because  it  was  identified  as 
the  best  operating  condition  in  [3].  The  experimentally  observed  test  time  was  approximately 
300  fis  before  large-amplitude  noise  contaminated  the  test  flow. 

Each  case  was  run  as  a  sequence  of  two  hour  jobs  on  the  Cray-YMP  with  the  flow  data 
being  written  every  0.5  ms  (simulation  time)  for  the  shock- compression  process  and  every 
0.1  ms  for  the  test-gas  expansion.  Table  2  shows  the  required  computer  resources  together 
with  the  name  of  the  mass-store  directory  where  the  data  files  have  been  archived. 

Two  calculations  were  performed.  Case  la  used  a  “coarse”  grid  of  1604  cells  axially  by 
20  cells  radially  while  case  lb  used  a  “fine”  grid  of  2406  x  30  cells.  (Refer  to  Appendix  A 
for  the  coarse  grid  results.)  The  computed  results,  including  the  timing  of  events  such  as 
the  passage  of  the  shock  and  contact  surfaces,  were  largely  the  same.  There  were,  however, 
some  differences  in  the  details  of  the  contact  surfaces  and  the  noise  introduced  to  the  test 
flow. 

Figure  5  summarizes  the  calculation  for  case  lb  as  a  wave  diagram.  The  data  points 
plotted  were  obtained  from  the  individual  plots  of  the  density  field  shown  later.  Note  the 
small  size  of  the  slug  of  compressed  acceleration  gas  (state  20).  The  position  of  the  long 
dashed  line  denoting  the  -downstream-tail  of  the  unsteady  expansion  of  the  test  gas  was 
estimated  by  its  arrival  at  x  =  24  m  while  the  short-dashed  line  represents  the  propagation 
of  the  leading  edge  of  the  reflected  expansion  and  the  asterisk  denotes  the  arrival  of  large- 
amplitude  noise  at  the  end  of  the  acceleration  tube  (x  ~  24  m).  This  reflected  expansion  is 
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the  result  of  the  interaction  between  the  upstream-head  of  the  unsteady  expansion  and  the 
test-gas/driver-gas  interface.  The  test  time  (state  5)  for  this  case  starts  after  the  passage  of 
the  test-gas/acceleration-gas  interface  and  stops  with  the  arrival  of  the  downstream  tail  of 
unsteady  expansion. 

3.2.1  Shock  Compression  of  the  Test  Gas 

Density  contours  for  the  shock  compression  of  the  test  gas  are  shown  in  Fig.  6.  For  each 
frame,  81  contours  have  been  plotted  over  the  range  —6.0  <  log^o  P  <  2.0.  Although  the 
shock  has  remained  essentially  planar,  the  driver-gas/test-gas  interface  has  become  quite 
distorted  hy  t  =  2.0  ms.  Boundary  layer  effects  are  present  (as  can  be  seen  by  the  contact 
surface  being  left  behind  along  the  tube  wall)  but  viscous  effects  are  not  well  resolved  in  this 
part  of  the  calculation. 

Of  special  interest  is  the  train  of  transverse  waves  following  the  interface.  These  transverse 
waves  show  up  clearly  in  the  contours  of  the  radial  velocity  as  shown  in  Fig.  7  where  it 
appears  that  the  disturbances  are  introduced  at  the  primary  diaphragm  station.  Note  that 
the  axial  scale  is  highly  compressed  and  the  actual  Mach  angles  are  much  shallower  than 
those  that  appear  in  the  figure.  Significant  noise  is  present  in  the  shocked  test  gas. 

Figure  8  shows  the  axial  variation  of  flow  properties  at  f  =  1.0  ms  together  with  the 
estimated  properties  for  a  purely  one-dimensional  situation.  Except  for  some  noise,  the 
levels  are  in  general  agreement.  Although  the  axisymmetric  simulation  exhibits  a  higher 
shock  speed  than  the  one-dimensional  case,  tables  3  and  4  show  that  the  computed  shock 
speed  is  still  lower  than  that  measured  in  the  experiment.  This  effect  was  discussed  in  [3] 
and  may  be  partly  due  to  a  low  estimate  for  the  initial  driver-gas  temperature.  Also,  note 
the  variation  in  shock  speed  over  this  phase  of  the  facility  operation  (see  table  4). 

3.2.2  Unsteady  Expansion  of  the  Test  Gas 

The  unsteady  expansion  of  the  test  gas  is  shown  in  Fig.  9.  Viscous  effects  are  more  important 
in  this  process  and  a  very  thick  boundary  layer  is  established  behind  the  shock  and  test- 
gas/acceleration-gas  interface.  The  separation  distance  of  the  shock  and  contact  remains 
essentially  constant  from  t  =  2.5  ms  to  t  =  4.0  ms  because  a  balance  has  been  reached 
between  the  mass  flow  of  acceleration  gas  entering  through  the  shock  and  the  mass  flow 
into  the  boundary  layer  [10].  This  hot  acceleration  gas  tends  to  accumulate  just  upstream 
of  the  test-gas/acceleration-gas  interface  and  form  a  large  bulge  in  the  boundary  layer  (at 
X  ~  23  m,<  =  4.5  ms).  It  is  not  until  the  last  frame  (t  =  5.0  ms)  that  any  significant  level 
of  noise  is  evident  in  the  expanding  test  gas. 
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Figure  10  shows  the  axial  variation  of  flow  properties  &t  t  =  4.4  ms.  There  are  very 
obvious  disturbances  in  the  region  of  the  expanded  driver  gas  (state  3)  and,  on  close  inspec¬ 
tion  of  the  pressure  distribution,  some  small-amplitude  disturbances  can  be  seen  near  the 
downstream  end  of  the  unsteady  expansion  of  the  test  gas  (at  a;  ~  22  m).  The  very  high 
temperature  of  the  shocked  acceleration  gas  is  interesting. 

Figures  11,  12  and  13  show  details  of  the  shock  and  test-gas/acceleration-gas  interface 
region.  The  bulge  seen  very  clearly  at  a:  ~  23  m,  t  =  4.5  ms  in  Fig.  9  is  no  longer  as 
obvious  in  these  “true-shape”  plots.  Plotting  the  velocity  field  in  a  frame  of  reference  which 
is  stationary  with  respect  to  the  shock  and  contact  surface  clearly  shows  the  bleeding  of 
shocked  acceleration  gas  out  of  region  20  and  into  the  boundary  layer.  Estimates  of  the 
separation  distance  between  the  shock  and  contact  surface  (i.e.  the  length  of  region  20)  can 
be  made  using  the  theory  presented  in  [10].  Using  the  expression 

(11) 

4p  \Pw,0  J  ‘^tu  'Uefl  Pw,0 
together  with  the  parameter  estimates 

/3  =  2.12,  =  C/a,io  =  5900  m/s,  Ue.o  =  1200  m/s. 


pe,o  =  1.0  X  10"'^  Icg/m^,  p^^  =  2.484  x  10  kg/m^^  =  2.012  x  10  ®  Pa.s, 

a  separation  length  of  0.2  m  is  obtained.  This  is  in  reasonable  agreement  with  the  length 
observed  in  Fig.  11  but  it  should  be  noted  that  the  value  of  L20  obtained  from  eqn.  (11)  is 
very  sensitive  to  the  estimated  parameters. 


3.2.3  Histories  of  the  Test  Flow  Properties 

Histories  of  the  centreline  test-flow  properties  near  the  exit  of  the  acceleration  tube  are 
shown  in  Fig.  14.  Features  that  can  be  identified  include  the  passage  of  the  shock  and  the 
test-gas/acceleration- gas  interface  at  times  t  =  4.45  ms  and  t  =  4.83  ms  respectively.  These 
events  are  followed  by  a  gradual  variation  in  properties  possibly  due  to  a  constriction  of  the 
effective  tube  diameter  associated  with  the  passage  of  the  bulge  (identified  in  Fig.  9).  The 
downstream  tail  of  the  unsteady  expansion  arrives  at  <  =  4.63  ms  and  average  pressure  and 
density  begin  to  rise  while  Mach  number  and  velocity  decline.  It  is  in  the  period  4.6  <  t  <  5.0 
that  a  “low-amplitude”  fluctuation  is  evident. 

The  fluctuations  increase  in  both  the  amplitude  and  the  frequency  with  the  arrival  of 
the  reflected  head  of  the  unsteady  expansion  at  t  =  5.04  ms.  This  wave  originates  with  the 
rupture  of  the  secondary  diaphragm  and  travels  along  a  1*2  —  02  characteristic  until  it  strikes 
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the  driver-gas/test-gas  interface.  There  is  both  transmission  into  the  driver  gas  and  partial 
reflection  downstream  along  &  u  +  a  characteristic  thus  admitting  more  disturbances  into 
the  expanding  test  gas.  The  arrival  time  for  this  wave  (at  the  exit  plane  of  the  acceleration 
tube)  may  be  estimated  as  described  in  [6].  We  start  at  x'  =  0,  t'  =  0  with  the  rupture 
of  the  secondary  diaphragm.  The  upstream  head  of  the  unsteady  expansion  intersects  the 
driver-gas/test-gas  interface  at  t'  =  Ati,  x'  =  (u2  —  0,2)  Ati.  Refer  to  tables  5  and  6  for 
particular  values.  The  propagation  time  of  the  wave  through  the  unsteady  expansion  may 
be  obtained  by  integration  along  a,  u  +  a  characteristic  for  which  the  Riemann  invariant  is 
[11] 


Us  + 


7  -  1 

Within  the  centred  expansion  we  have 


Qs  —  U  —  U  + 


7-1 


X 

—  =  u  —  a 
f 


(12) 


(13) 


and  this  may  be  substituted  into  equation  (12)  to  obtain 


a 

u 


7  -  1 
7  +  1 


2  x'\ 

7-HV 


(14) 


The  integration  may  be  terminated  at  the  end  of  the  expansion  (a  =  C5)  or  at  the  end  of  the 
acceleration  tube,  whichever  is  sooner.  If  the  reflected  wave  reaches  the  downstream  end  of 
the  expansion  before  it  reaches  the  end  of  the  tube,  it  continues  to  propagate  into  the  test 
flow  region  (20)  with  speed  115  +  05.  The  time  of  arrival  as  computed  with  this  approach 
closely  matches  the  onset  of  large-amplitude  noise  in  the  exit  flow  histories. 


A  comparison  of  computed  test  flow  properties  with  the  measurements  of  Shinn  and 
Miller  [3]  is  shown  in  Fig.  15.  The  overall  level  of  the  computed  wall  pressure  is  a  little  low, 
possibly  because  of  the  low  estimate  for  the  initial  temperature  in  the  driver  tube.  Also,  the 
computed  wall  pressure  is  smoother  than  the  experimental  wall  pressure  but  this  may  be 
partly  due  to  numerical  dissipation  and  partly  due  to  noisy  measurement  techniques  in  the 
experiment.  The  good  news  is  that  arrival  of  the  downstream  tail  of  the  unsteady  expansion 
(0.2  ms  after  the  shock)  agrees  well. 


The  other  frames  in  Fig.  15  show  pitot  pressures.  The  computed  pitot  pressure  at  the 
centreline  is  again  a  little  lower  than  the  experimental  value  and  is  also  less  noisy  for  the 
first  0.4  ms  after  shock  arrival.  However,  there  is  good  agreement  in  that  there  appear  to 
be  two  types  of  noise:  a  small-amplitude  noise  before  the  arrival  of  the  reflected  head  of  the 
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unsteady  expansion  and  a  large-amplitude  noise  after  the  passage  of  the  reflected  expansion. 
This  is  one  of  the  main  results  of  this  study. 

Other  features  evident  in  the  pitot  histories  include  the  large  dip  at  <  =  5.0  ms,  r  = 
52.5  mm  and  the  consistently  low  values  at  r  =  71  mm.  The  large  dip  is  associated  with 
the  passage  of  the  bulge  in  the  boundary  layer  while  the  r  =  71  mm  position  remains  within 
the  boundary  layer  after  the  passage  of  the  test-gas/acceleration-gas  interface. 

3.3  Case  2:  Pi  =  20.62  kPa,  Piq  —  16  Pa 

Having  successfully  simulated  the  arrival  of  noise  in  the  test  flow  at  Pi  =  3.45  kPa,  we 
decided  to  try  a  higher  pressure  of  Pi  =  20.62  kPa  where  the  experimental  measurements 
indicated  that  large-amplitude  disturbances  appear  very  early  in  the  test  flow.  Fig.  16  shows 
the  wave  diagram  for  this  case.  Note  that  the  downstream  tail  of  the  unsteady  expansion 
arrives  at  the  tube  exit  shortly  after  the  test  flow  begins  and  that  the  reflected  head  of  the 
unsteady  expansion  is  expected  to  arrive  much  later. 

The  shock  compression  of  the  test  gas  is  shown  in  Fig.  17.  Because  of  the  higher  initial 
intermediate-tube  pressure,  the  shock  has  a  lower  speed  than  in  case  1.  Transverse  waves 
are  still  present  in  the  expanded  driver  gas  but  they  seem  to  be  weaker.  Also,  there  is  some 
noise  in  the  compressed  test  gas. 

The  unsteady  expansion  of  the  test  gas  is  shown  if  Fig.  18.  The  process  appears  similar 
to  that  in  case  1  except  that  there  is  very  little  noise  evident  in  the  expanded  test  gas  (even 
at  f  =  5.0  ms).  This  lack  of  noise  is  confirmed  in  Fig.  19  and  Fig.  20.  In  Fig.  20,  the  arrival 
of  the  downstream-tail  of  the  unsteady  expansion  can  be  clearly  seen  at  f  ~  4.8  ms  while 
the  reflected  head  of  the  unsteady  expansion  is  not  expected  to  arrive  at  the  tube  exit  until 
t  ~  5.84  ms  (see  table  6). 

Fig.  21  compares  the  computed  pressures  at  the  exit  of  the  acceleratior  tube  with 
experimental  measurements.  The  lack  of  noise  in  the  laminar  solution  suggests  that  the 
disturbances  in  the  experimental  test  flow  are  caused  by  boundary  layer  transition  or  turbu¬ 
lence  (which  is  not  modelled  in  these  simulations)  rather  than  propagation  of  noise  through 
the  driver-gas/ test-gas  interface. 

3.4  Case  3:  Pi  =  1.74  kPa,  Piq  =  16  Pa 

Given  the  failure  of  the  Pi  =  20.62  kPa  case  to  exhibit  large  disturbances  in  the  simulated 
test  flow,  we  decided  to  try  a  case  with  lower  intermediate-tube  pressure  than  case  1.  Figures 
22  -  27  show  the  results  of  a  simulation  with  Pi  =  1.74  kPa,  Pio  =  16  Pa. 
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Table  4  shows  that  the  shock  speed  during  the  initial  compression  of  the  test  gas  is  slightly 
higher  than  for  case  1  but  that  the  final  shock  speed  during  the  unsteady  expansion  process 
is  slightly  lower.  The  shock  compression  of  the  test  gas  (Fig.  23)  appears  to  be  noisier  than 
for  case  1  and  the  centreline  histories  (Fig.  26)  indicate  that  small-amplitude  disturbances 
appear  very  early  in  the  test  flow.  This  seems  to  be  consistent  with  the  experim^'  .^al  pitot 
pressure  results  shown  in  Fig.  27.  The  arrival  of  the  reflected  head  of  the  unsteady  expansion 
(t  ~  5.1  ms)  can  be  seen  as  a  decrease  in  the  growth  rate  of  pitot  pressure  and,  in  the 
simulation,  is  accompanied  by  larger-amplitude  and  higher-frequency  disturbances  which 
are  not  obvious  in  the  experimental  measurements. 

3.5  Case  4:  Pi  =  3.45  kPa,  Piq  =  52.6  Pa 

For  large  values  of  acceleration-tube  filling  pressures,  Shinn  and  Miller  [3]  also  observed  a 
large  dip  in  the  pitot  pressure  shortly  after  the  start  of  the  test  flow.  The  phenomenon  was 
attributed  to  boundary  layer  transition.  We  decided  to  try  simulating  such  an  operating 
condition  to  see  if  the  bulge  in  the  boundary  layer  (following  the  test-gas/acceleration-gas 
contact  surface)  was  a  contributing  factor. 

The  results  for  the  laminar-flow  simulation  with  Pi  =  3.45  kPa,  Pio  =  52.6  Pa  are 
shown  in  figures  28  -  32.  The  shock  compression  phase  is  the  same  as  for  case  la  but, 
with  the  larger  mass  of  acceleration  gas,  the  separation  distance  between  the  shock  and  the 
contact  surface  is  increased  and  equilibrium  is  attained  hy  t  =  2.8  ms.  The  histories  of 
the  centreline  properties  shown  in  Fig.  31  indicate  that  small-amplitude  noise  is  present 
in  both  the  acceleration  gas  and  test  gas  up  until  the  arrival  of  the  reflected  head  of  the 
unsteady  expansion  at  ^  ~  4.9  ms.  Although  the  large-amplitude  noise  after  the  passage  of 
the  reflected  expansion  seems  to  be  consistent  with  the  experimental  measurements  shown 
in  Fig.  32,  there  is  no  evidence  of  the  large  pitot-pressure  dip  in  the  laminar  simulation. 
This  is  consistent  with  the  hypothesis  [3]  that  the  dip  is  caused  by  boundary  layer  transition. 
Note  that  this  case  was  computed  on  a  coarse  mesh  which  may  result  in  the  attenuation  of 
the  high-frequency  components  of  the  noise  (see  appendix  A). 

4  Concluding  Remarks 

From  the  computational  point  of  view,  we  can  now  make  reliable  simulations  for  laminar 
aixisymmetric  flow  of  a  perfect  gas  in  an  expansion  tube.  Viscous  effects  were  adequately  re¬ 
solved  in  the  acceleration  tube  but  they  were  not  well  resolved  in  either  the  intermediate  tube 
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or  the  driver  tube.  This  is  not  considered  to  be  a  major  problem  as  inviscid  multidimensional 
flow  features  seem  to  dominate  the  performance  of  these  sections  of  the  facility. 

The  numerical  technique  is  capable  of  capturing  strong  shocks  and  contact  surfaces 
with  large  temperature  jumps,  however  the  contact  surfaces  tend  to  be  more  (numerically) 
smeared  them  the  shocks.  The  net  result  is  that  the  estimates  of  wave  propagation  times  and 
test  flow  properties  appear  to  be  in  reasonable  agreement  with  experimental  measurements. 
However,  the  issue  of  reliably  capturing  the  high-frequency  noise  is  not  fully  settled  and  fur¬ 
ther  computational  experiments  using  alternative  higher-order  schemes  (such  as  ENO  [12]) 
are  required. 

We  deliberately  selected  a  laminar  perfect-gas  flow  to  keep  this  study  manageable  but 
there  is  an  immediate  need  to  perform  simulations  with  nonideal  thermal  and  chemical  effects 
with  turbulence  in  three-dimensions.  Diaphragm  dynamics  is  another  important  item  which 
needs  to  be  considered  for  future  study.  In  the  near  future,  it  is  unlikely  that  we  could 
include  all  of  these  features  but  the  modelling  of  turbulence  and  the  diaphragm  dynamics 
are  key  items  required  to  simulate  the  expansion  tube  over  a  larger  range  of  useful  operating 
conditions. 

With  respect  to  facility  design  and  operation,  the  simulations  show  some  important 
features.  Firstly,  the  discontinuity  at  the  primary  diaphragm  station  is  sufficient  to  produce 
disturbances  which  may  then  propagate  into  the  test  flow.  Secondly,  the  simulations  identify 
two  categories  of  noise  in  the  test  flow.  One  is  a  small-amplitude  disturbance  which  may  be 
seen  before  the  arrival  of  the  reflected  head  of  the  unsteady  expansion  while  the  other  is  a 
large-amplitude  noise  associated  with  the  passage  of  the  reflected  expansion. 

This  distinction  explains  the  experimental  observation  [4,  5]  that  increasing  the  acceleration- 
tube  pressure  Piq  while  keeping  the  intermediate-tube  pressure  Pi  constant  results  in  higher 
frequency  noise  in  the  test  flow.  The  increase  in  Pio  results  in  a  lower  shock  speed  in  the 
acceleration  tube,  while  the  fixed  Pi  means  that  the  arrival  time  of  the  reflected  expansion 
will  be  essentially  constant.  Thus,  at  high  values  of  Fio,  we  expect  the  reflected  expansion  to 
arrive  very  soon  after  the  shock  and  introduce  high-frequency  disturbances  to  the  test  flow. 

To  avoid  the  large-amplitude  noise,  it  should  be  relatively  simple  to  select  operating 
conditions  or  facility  configurations  which  delay  the  arrival  of  the  reflected  expansion.  For 
example,  the  intermediate  tube  could  be  lengthened  while  fixing  the  initial  gas  pressures. 

It  is  not  immediately  obvious  how  to  avoid  the  small-amplitude  disturbances  observed 
before  the  arrival  of  the  unsteady  expansion.  As  can  be  seen  in  Fig.  33  these  disturbances  do 
not  appear  in  the  test  flow  for  all  operating  conditions.  It  was  hoped  that  the  acoustic  theory 
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of  Pauli  would  provide  some  guide  to  the  correct  selection  of  good  operating  conditions.  The 
basic  result  of  that  theory  is  that  noise  crossing  the  driver- gas/ test-gas  interface  is  attenuated 
if  the  speed  of  sound  in  the  test  gas  is  much  higher  than  the  speed  of  sound  in  the  driver  gas. 
This  jump  in  sound  speeds  is  present  for  each  of  the  cases  studied  here  but,  in  contrast  to 
Pauli’s  theory,  the  case  with  the  least  noise  also  had  the  smallest  jump  in  sound  speed  across 
the  driver-gas/test-gas  interface  (see  table  5),  We  suggest  that  a  number  of  other  operating 
conditions  be  simulated  before  drawing  any  firm  conclusions. 
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Appendix  A 


Coarse  Grid  Calculation,  Case  la. 

In  order  to  partially  address  the  issue  of  grid  sensitivity,  we  simulated  the  Pi  =  3.45  kPa, 
Pio  =  16  Pa  condition  on  a  coarse  (1604  x  20)  grid.  This  case  is  labelled  la  in  tables  1  - 
3.  The  results  are  shown  in  figures  34  -  37.  Except  for  fine-scale  detail,  the  results  are 
essentially  the  same  as  for  the  2406  x  30  (fine)  grid  simulation.  The  most  notable  difference 
is  shown  in  Fig.  37  which  compares  the  test  flow  pressures  at  the  exit  of  tne  acceleration 
tube.  The  coarse  grid  has  attenuated  the  small-amplitude  noise  seen  before  the  arrival  of 
the  reflected  head  of  the  unsteady  expansion  and  has  filtered  out  most  of  the  high-frequency 
component  of  the  noise  after  the  passage  of  the  reflected  expansion  (at  t  ~  5.0  ms).  The 
(relative)  timing  of  events  and  the  mean  levels  of  the  pressures  are  unchanged  but  there  is  a 
slight  time  shift  in  the  plots  in  order  to  offset  a  timing  mistake  in  resetting  the  acceleration 
tube  initial  conditions. 
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Case 

Pi 

(kPa) 

{kglm^) 

Pio 

{Pa) 

p\o 

(%/m3) 

la 

3.45 

5.537  X  10-3 

16 

2.568  X  10-® 

lb 

3.45 

5.537  X  10-3 

16 

2.568  X  10-® 

2 

20.62 

3.309  X  10-2 

16 

2.568  X  10-® 

3 

1.74 

2.792  X  10-3 

16 

2.568  X  10-3 

4 

3.45 

5.537  X  10-3 

52.6 

8.442  X  10-® 

Table  1:  Initial  States. 


Case 

CFL 

cpu  —  time 
(hrs) 

memory 
(M  words) 

^  final 

{ms) 

grid 

mass-store 

directory 

la 

0.5 

12.5 

3.2 

5.4 

1604  X  20 

exp-D2 

lb 

0.5-0.6 

42 

6.8 

5.4 

2406  X  30 

exp-B 

2 

0.5-0.6 

41.6 

6.8 

5.3 

2406  X  30 

exp-C 

3 

0.5 

49 

6.8 

5.4 

2406  X  30 

exp-F 

4 

0.5 

6.5 

3.2 

5.4 

1604  X  20 

exp-G 

Table  2:  Computational  resources. 


Table  3:  Experimentally  measured  shock  speeds  and  pressures  from  Shinn  and  Miller  (1978). 


Table  4:  Computed  shock  speeds  and  pressures. 
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Table  5:  Flow  properties  across  the  unsteady  expansion  of  the  test  gas. 


Case 

^rupture 

(ms) 

Ati 

(ms) 

Atpropagate 

(ms) 

^arrival 

{ms) 

lb 

1.00 

2.96 

5.04 

a 

1.15 

3.44 

5.84 

n 

2.00 

0.95 

2.82 

4.82 

2.08 

0.95 

2.95 

5.03 

Table  6:  Event  times  associated  with  the  reflection  of  the  unsteady  expansion. 
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Primary  diaphragn  Secondary 

(single  or  double)  diaphraon 

Unsteady 


Figure  1:  Conceptual  wave  diagram  for  the  operation  of  the  NASA  Langley  expansion  tube. 
The  numbers  refer  to  the  states  of  the  gases  as  defined  by  Trimpi  (1962). 
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Figure  3;  Comparison  of  the  density  fields  at  <  =  2.0  ms  for  case  2:  (a)  1604  x  20  grid,  Euler 
time-stepping;  (b)  1604  x  20  grid,  predictor-corrector  time-stepping;  (c)  2406  x  30  grid,  Euler 
time-stepping. 


Figure  4:  Expansion  tube  geometry  as  modelled  in  the  simulations.  Dimensions  are  in 
metres. 


secondary  diaphragm  rupture 


Figure  5:  Wave  diagram  for  case  1.  Q  =  shock;  A  =  driver-gas/test-gas  interface;  +  = 
test-gas/acceleration- gas  interface;  long  dash  =  downstream  tail  of  the  unsteady  expansion; 
o  =  upstream  head  of  the  unsteady  expansion;  short  deish  =  reflected  expansion. 


Pi  =  3.45  kPa,  Pm  =  16  Pa 


primary  diaphragm 


t  =  0.5  ms 


contact— 1  /-shock  .  i  n 

\  jf  t  =  1.0  ms 


1.5  ms 


Figure  6:  Density  contours  showing  the  shock  compression  of  the  test  gas  for  case  1.  The 
flow  states  are  labelled  as  per  Fig.  1, 
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Pi  =  3.45  kPa^  Pio  =  16  Pa 


Expansion  Tube  (B),  t  =  1.0ms,  iy  =  1 


X,  m 


X,  m 


Figure  8:  Flow  properties  along  the  axis  during  shock  compression  (case  1,  <  =  1.0  ms).  The 
solid  line  denotes  the  ideal  one-dimensional  solution. 
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P\  —  3.45  kPa,  P^^y  IG  Pa 


Pi  =  3.45  kPa,  Pio  =  16  Pa 


(b) 


Figure  10:  Properties  along  the  axis  at  f  =  4.4  ms  for  case  1:  (a)  density  contours;  (b)  other 
properties  as  labelled. 
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P\  =  3.45  kPa,  Pio  =  16  Pa 


Expansion  Tube  (B),  t  =  4.4ms.  iy  =  1 


X,  m  X,  m 


22.0  22.4  22.8  23.2  23.6  24.0  22.0  22.4  22.8  23.2  23.6  24.0 

X,  m  X,  m 


Figure  11:  Properties  along  the  axis  and  near  the  shock  at  t  =  4.4  ms,  case  1. 
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Pi  =  3.45  kPa,  Pio  =  16  Pa 


Figure  12:  True-shape  plots  of  the  flow  near  the  shock  at  t  =  4.4  ms,  case  1:  (a)  density 
contours;  (b)  pressure  contours;  (c)  specific-energy  contours. 


T 


Pi  =  3.45  kPa,  Pio  =  16  Pa 


X,  m 


Figure  16:  Wave  diagram  for  case  2.  Q  =  shock;  A  =  driver- gas/test-gas  interface;  -f  = 
test-gM/ acceleration-gas  interface;  long  dash  =  downstream  tail  of  the  unsteady  expansion; 
o  =  upstream  head  of  the  unsteady  expansion;  short  dash  =  reflected  expansion. 


Figure  17:  Density  contours  showing  the  shock  compression  of  the  test  gas  for  case  2 
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Pi  =  20.62  kPa,  Pio  =  16  Pa 


CFD  Results 


Experimental  Results 


Figure  21:  Comparison  of  wall  pressure  and  pitot  pressure  with  experiment  (case  2,  x 
24  m).  The  experimental  data  is  taken  from  Shinn  and  Miller  (1978). 
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X,  m 


Figure  22:  Wave  diagram  for  case  3.  0  =  shock;  A  =  driver-gas/test-gas  interface;  +  = 
test-gas/acceleration-gas  interface;  long  dash  =  downstream  tail  of  the  unsteady  expansion; 
o  =  upstream  head  of  the  unsteady  expansion;  short  dash  =  reflected  expansion. 


Pi  =  1.74  kPa,  Pio  =  16  Pa 


0.10 


Pi  =  1.74  kPa,  Pio  =  16  Pa 


Figure  25:  Properties  along  the  axis  at  <  =  4.4  ma  for  case  3:  (a)  density  contours;  (b)  other 
properties  eis  labelled. 
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Figure  27:  Comparison  of  wall  pressure  and  pitot  pressure  with  experiment  (case 
24  m).  The  experimental  data  is  taken  from  Shinn  and  Miller  (1978). 


X,  m 


Figure  28:  Wave  diagram  for  case  4.  Q  =  shock;  A  =  driver-gas/test-gas  interface;  +  = 
test-gas/acceleration-gas  interface;  long  dash  =  downstream  tail  of  the  unsteady  expansion 
0  =  upstream  head  of  the  unsteady  expansion;  short  dash  =  reflected  expansion. 
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Density  contours  showing  the  unsteady  expansion  of  the  test  gas  for  case  4. 
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=  3.45  kPa,  Pio  =  52.6  Pa 


0.10  f- 


0.00 


la) 


(b) 


Figure  30:  Properties  along  the  axis  at  <  =  4.5  ms  for  case  4:  (a)  density  contours;  (b)  other 
properties  as  labelled. 


Pi  =  3.45  kPa,  Pio  =  52.6  Pa 
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CFD  Results 


Experimental  Results 


4.4 


4.8 


5.2 


100 


4  T 


kPa  2  I 


1  • 

0  - • — 

0  0.2 

t,ms 


■H 

0.4 


Figure  32:  Comparison  of  wall  pressure  and  pitot  pressure  with  experiment  (case  4,  x 
24  m).  The  experimental  data  is  taken  from  Shinn  and  Miller  (1978). 


.9 


2000 


Case  lb,  t  =  4.4  ms 


Figure  33:  Axial  variation  of  the  sound  speed  f< 
at  the  acceleration-tube  exit. 
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0  5  10  15  ^  ^  20  25 


Figure  35.  Density  contours  showing  the  unsteady  expansion  of  the  test  gas  for  case  la. 
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0  5  10  15  20  25 


X,  m 


Figure  36:  Properties  along  the  axis  at  ^  =  4.4  ms  for  case  la:  (a)  density  contours;  (b) 
other  properties  as  labelled. 
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